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ABSTRACT 

Single-molecule techniques facilitate analysis of 
mechanical transitions within nucleic acids and 
proteins. Here, we describe an integrated fluores- 
cence and magnetic tweezers instrument that 
permits detection of nanometer-scale DNA struc- 
tural rearrangements together with the application 
of a wide range of stretching forces to individual 
DNA molecules. We have analyzed the force- 
dependent equilibrium and rate constants for 
telomere DNA G-quadruplex (GQ) folding and un- 
folding, and have determined the location of the 
transition state barrier along the well-defined 
DNA-stretching reaction coordinate. Our results 
reveal the mechanical unfolding pathway of the 
telomere DNA GQ is characterized by a short 
distance (<1 nm) to the transition state for the 
unfolding reaction. This mechanical unfolding 
response reflects a critical contribution of long- 
range interactions to the global stability of the 
GQ fold, and suggests that telomere-associated 
proteins need only disrupt a few base pairs to 
destabilize GQ structures. Comparison of the GQ 
unfolded state with a single-stranded polyT DNA 
revealed the unfolded GQ exhibits a compacted 
non-native conformation reminiscent of the protein 
molten globule. We expect the capacity to interro- 
gate macromolecular structural transitions with 
high spatial resolution under conditions of low 
forces will have broad application in analyses of 
nucleic acid and protein folding. 

INTRODUCTION 

Model RNA and DNA hairpins have been extensively 
characterized by optical trapping methods, providing a 



detailed view of the folding energy landscapes of these 
fundamental nucleic acid structures (1-5). However, the 
ability to achieve high spatial resolution using force spec- 
troscopy relies on the application of relatively large 
stretching forces (>10pN) to suppress the measurement 
noise introduced by the long flexible DNA handles used 
to attach the structure of interest to micron-scale beads 
held in the optical trap. Recently, a study of the force- 
dependent structural dynamics of single Holliday junc- 
tions was reported, which paired an optical trap to 
apply precisely calibrated stretching forces with single- 
molecule Fdrster resonance energy transfer (smFRET) 
to monitor DNA structural dynamics (6). This fluores- 
cence-force method provided a powerful tool for 
probing sub-nanometer-scale structural rearrangements 
within single Holliday junctions at very low stretching 
forces (<lpN) over short periods. In addition, several 
groups have reported measurements that combine 
smFRET with a magnetic tweezers apparatus (7,8). The 
use of magnets to apply mechanical loads to individual 
DNA molecules has several potential advantages over 
optical traps. For example, the combination of a high- 
power optical trapping laser with single-molecule fluores- 
cence is technically difficult to implement, typically 
requiring interlacing of the FRET excitation and optical 
trapping beams to avoid rapid photo damage of the FRET 
probes by the high-power trapping laser (9,10). Moreover, 
the sensitivity of optical traps to mechanical drift makes 
the application of low stretching forces (<lpN) over 
extended periods far more challenging than with a 
simple magnetic tweezers system. Here we describe an 
integrated fluorescence and magnetic tweezers microscope 
capable of measuring nanometer-scale structural transi- 
tions in single DNA molecules at low stretching forces. 
We demonstrate the utility of this approach by analyzing 
the mechanical unfolding pathway of a model human 
telomere DNA substrate. 

Telomeres are specialized chromatin structures that 
protect linear ends of eukaryotic chromosomes from 
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aberrant DNA processing by DNA damage repair ma- 
chinery (11). The foundation of human telomere structure 
is a long stretch of double-stranded DNA composed of a 
hexanucleotide DNA repeat sequence (T 2 AG 3 ). In 
addition, all telomeres terminate with a 3' single-stranded 
G-rich DNA tail, which has the capacity to fold into a 
unique secondary structure called a G-quadruplex (GQ). 
Human telomere DNA GQs are proposed to play a 
central role in telomere homeostasis, and small-molecule 
ligands that selectively bind and stabilize telomere DNA 
GQs have shown promise as anti-cancer drugs (12,13). 
Thus, intensive efforts have been made to better under- 
stand the structure and function of telomere DNA GQs. 
The first solution structure of a human telomere GQ 
revaled a fundamental structural architecture in which 
guanine bases are hydrogen bonded in a planar quartet 
geometry and may coordinate a single centrally located 
monovalent metal ion (Figure 1A, top left) (14). Three 
adjacent intra-molecular G-quartets may interact via 
stacking interactions and are topologically linked by 
short intervening DNA loop sequences (Figure 1A, 
bottom left). Moreover, the folding properties of 
telomere DNA GQs vary with the presence of different 
monovalent cations. Na + ions predominantly promote 
the formation of an anti-parallel GQ conformation (14), 
whereas GQ DNA crystals formed in the presence of K + 
ions revealed a distinct parallel GQ folding topology (15). 
More recent solution studies have demonstrated that 
multiple GQ topologies coexist in the presence of K + , 
including the anti-parallel, parallel and several hybrid 
forms (16-19). 

The structural and dynamic properties of telomere 
DNA GQs have been studied using smFRET (20-22). 
These experiments revealed that a particular GQ topo- 
logical fold must transit through an obligatory unfolded 
intermediate to isomerize into a distinct GQ fold. More 
recently, the rupture force distribution of single telomere 
DNA GQs has been analyzed using atomic force micros- 
copy and optical trapping, providing a direct measure- 
ment of telomere DNA GQ mechanical stability (23-26); 
however, these force spectroscopy studies did not analyze 
the force dependence of telomere DNA GQ folding/un- 
folding at equilibrium. Here, we report a detailed analysis 
of the force dependence of telomere DNA GQ folding and 
unfolding at equilibrium under Na + folding conditions. 
Our results demonstrate that the structural equilibrium 
between the unfolded and Na + -induced telomere DNA 
GQ structure is highly sensitive to forces between ~1 
and 8pN. Analysis of the force-dependent rate constants 
for folding and unfolding provides a direct measurement 
of the position of the transition state barrier for telomere 
DNA GQ unfolding along the DNA-stretching reac- 
tion coordinate. Interestingly, in contrast to other DNA 
secondary structures characterized by force spectroscopy, 
we find telomere DNA GQs exhibit a very short distance 
(<1 nm) to the transition state barrier for unfolding. This 
unfolding behavior indicates telomere DNA GQ structure 
is significantly stabilized by long-range contacts, and once 
these contacts are disrupted, the entire GQ fold readily 
dissolves. 



MATERIALS AND METHODS 

DNA oligonucleotides 

All DNA oligonucleotides were purchased from 
Integrated DNA Technologies, Inc. The sequences of 
all DNAs used in the study are listed in Supplementary 
Table SI. 

Dye labeling of DNA oligonucleotides 

The biotinylated DNA handle was labeled at the amino 
modification on C6 of T3 with monoreactive Cy5 (GE 
Healthcare), and the EcoRI DNA Handle was labeled at 
the amino modification on C6 of T16 with monoreactive 
Cy3 (GE Healthcare). Dye-labeled fragments were EtOH 
precipitated and purified by reverse-phase chromatog- 
raphy using a C8 column (Agilent, Eclipse XDB-C8) 
on an AKTA purifier. After HPLC purification samples 
were EtOH precipitated, and resuspended in ddH 2 0. 
DNA concentrations were determined using a 
NanoDrop 2000C (Thermo Scientific). 

DNA annealing reactions 

The biotinylated, Cy5-labeled and Cy3-labeled strands 
were annealed to the Tel24 abasic, 15R60/T8 hairpin 
abasic, or polyT abasic sequence by heating to 95°C for 
4 min, followed by slow cooling to room temperature 
(over several hours) in the presence of a buffer (50 mM 
Tris-HCl, pH 8) containing either 100 mM NaCl or 
lOOmM KC1. 

DNA molecule for integrated fluorescence and magnetic 
tweezers measurements 

DNA molecules were constructed by ligation of three pre- 
cursor DNA fragments: the annealed biotinylated 
smFRET fragment, a digoxigenin-modified DNA linker 
fragment and lambda DNA handle. The preparation of 
desired Cy5- and Cy3-labeled DNA insert is described 
earlier in the text. The digoxigenin-labeled DNA linker 
is synthesized by PCR reaction with a pUC19 template 
using primers flanking the multiple cloning site. PCR re- 
actions were set up using a dNTP mixture containing a 1 :4 
molar ratio of digoxigenin-1 1-dUTP (Roche): dTTP and 
digested with BamHI. The 15 721-bp lambda DNA handle 
was prepared by BamHI and EcoRI enzymatic digestion 
and agarose gel purification. Ligation reactions were 
set up with a 1:1:1 molar ratio of the desired Cy5- and 
Cy3-labeled DNA insert harboring an EcoRI sticky end, 
the purified EcoRI/BamHI lambda DNA handle and the 
BamHI digoxigenin-modified DNA linker fragment. 
Ligation reactions were run overnight at 16°C in the 
presence of T4 DNA ligase (NEB). The final ligation 
product was heat inactivated at 65°C to ensure the T4 
DNA ligase did not remain associated with the DNA. 

Single-molecule FRET measurements in the 
absence of force 

For these experiments, quartz slides (Finkenbeiner Inc.) 
were cleaned by sonicating for 20 min in 10% w/v 
Alconox, 20 min in ddH 2 0, 20 min in acetone, 20 min 
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in 1 M KOH and then 20 min in fresh 1 M KOH. Slides 
are then rinsed with ddH 2 0 and dried under nitrogen, 
followed by flame cleaning with a propane torch for ~2 
min. Sample chambers were prepared by sandwiching 
pieces of parafilm between the quartz slide and a 
plasma-cleaned cover glass (Harrik Plasma Cleaner), 
cover glasses from TED Pella, Inc. (Prod. No. 260146). 
Chambers were heated on a hot block (95°C) for 1 min to 
seal the parafilm to the glass. Channels (~10uL total 
volume) were treated with 35 uL of 1 mg/mL of 
biotinylated bovine serum albumin (Sigma Cat. 
No.A8549) for 5 min, washed with 100 uL of T50 buffer 
[10 mM Tris-HCl (pH 8), 50 mM NaCl], incubated with 
50 uL of 0.2mg/mL streptavidin (Invitrogen Cat. #S888) 
and then washed with 100 uL of T50. Channels were then 
equilibrated with buffer matching the desired experimental 
condition (see main text). 

Next, 100 uL of 5-10 pM annealed fluorescent and 
biotin-labeled smFRET DNA was deposited onto the 
streptavidin-coated quartz slide in the desired buffer con- 
dition. After 5-min incubation, desired buffer was flushed 
to the slide to remove the unbound DNA. Data were col- 
lected in imaging buffer containing 50 mM Tris (pH 8), 
100 mM KC1 or NaCl, 0.4% (w/v) D-glucose, O.lmg/ml 
glucose oxidase (Sigma Cat.# G2133-250KU), and 
0.02mg/ml catalase (CalBiochem Cat.#219001) and 
saturated with Trolox (Aldrich Cat.#23,881-3). Data 
were acquired using a green laser (532 nm, Laserglow, 
Inc.) and prism-type total internal reflection microscopy 
on an inverted Olympus 1X71 microscope equipped with 
an Andor IXON (897) CCD camera with 100-ms integra- 
tion time. 

Single-molecule FRET analysis 

Raw movie files were analyzed using in-house-written 
software available on request (IDL and MATLAB). 
FRET is defined as the efficiency of energy transfer 
between acceptor and donor dye, Ia/(Ia + Id)- Ia is the 
acceptor intensity, and I D is the donor intensity. FRET 
histograms were compiled by combining data from ~100 
individual single-molecule trajectories. The hidden 
Markov modeling was performed using the HaMMy 
(27) software program. First, dye intensity traces used in 
the HaMMy analysis were normalized after photo 
bleaching, and then the period before photo bleaching 
was cropped to be used in the analysis. Because individual 
traces did not typically possess a sufficient number of tran- 
sitions to produce a reliable HaMMy fit, we opted to 
stitch all the normalized and cropped trajectories into a 
single trace using in-house-written MATLAB software. 
This compiled trace was then fit using the HaMMy 
software. The idealized FRET trace that was produced 
by the HaMMy fitting was then parsed back into individ- 
ual traces using in-house MATLAB software to prevent 
the inclusion of artificial FRET transitions at the trace 
stitch points. The individual HaMMy fits were then used 
to generate a transition density plot using in-house- 
written MATLAB software. To fit the smFRET distribu- 
tions, data were binned (bin size was FRET = 0.03) and fit 
with multiple Gaussian functions using Origin software. 



For dwell time analysis, idealized traces produced by 
fitting the single-molecule trajectories using the HaMMy 
software package were used to calculate the average 
lifetime of the folded or unfolded states. The smFRET 
histograms for the Na + and K + titrations, as well as in 
the presence of Mg 2+ , were constructed by averaging the 
observed FRET for each molecule over 2 s (20 frames at 
100-ms integration time) and binning the data with a bin 
size of FRET = 0.03. 

Integrated FRET and magnetic tweezers measurements 

The DNA tethers for FRET-magnetic tweezers measure- 
ments were attached to anti-digoxigenin-coated magnetic 
beads. Anti-digoxigenin beads were made using EDC 
coupling chemistry (Pierce CAS# 22980) using 
carboxylated magnetic beads [Dynal, either 1-um 
(Invitrogen Cat.#650.11) or 2.8-um (Invitrogen 
Cat.#143.05D) diameter]. After incubating 10 min at 4°C, 
the DNA-magnetic bead mixture was washed to remove 
the excess DNA fragments before being deposited onto 
the streptavidin-coated cover glass (#1.5, Ted Pella) 
prepared using the same cleaning procedure described 
earlier in the text. After incubating 30 min at room tem- 
perature to allow the DNA-magnetic bead complexes to 
settle to the surface, 100 uL of 50 mM Tris with 100 mM 
Na + (for GQ experiments) or 200 mM K + (for control 
DNA hairpin experiments) was flowed onto the slide by 
gravity to remove the unattached DNA-linked magnetic 
beads. 

To calibrate the stretching forces being applied to each 
DNA tether, forces were measured across a wide range of 
positions of the magnet assembly held above the sample 
chamber. Force measurements were determined using the 
expression F = Lk B T/<x 2 >, where F is the force in 
pico-Newtons (pN), L is the DNA tether length at a par- 
ticular magnet height, k B T is thermal energy (4. 1 pN x nm) 
at room temperature and </ 2 > is the variance in the bead 
position in x-axis. The DNA tether length was measured by 
calibrating the diffraction ring pattern of the magnetic 
bead using a piezo-controlled objective positioning device 
(Mad City Labs) as previously described (28). Separate 
calibrations were made for the 1-um and 2.8-um 
magnetic beads. The force was found to decrease exponen- 
tially as the distance of the magnets was increased from the 
sample chamber as described (28). The magnet height 
versus force plots were fit with single exponential decay 
functions, and the parameters from these fits were then 
used to calculate the forces applied to subsequent experi- 
mental setups using the same-size bead. The bead manu- 
facturer estimates a deviation in the radius (r) of the beads 
on the order of 1-2%, and the force scales with the volume 
of the bead (r 3 ); thus, we estimate the systematic error of 
our force measurements based on an imperfect calibration 
to be on the order of 10%. 

Data were collected in imaging buffer containing 
100 mM Na + or 200 mM K + salt, 0.4% (w/v) D-glucose, 
O.lmg/ml glucose oxidase (Sigma Cat.# G2133-250KU), 
and 0.02mg/ml catalase (CalBiochem Cat.#2 19001) and 
saturated with Trolox (Aldrich Cat.#23,881-3). Data 
were acquired using objective-type total internal reflection 
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microscopy with a green laser (532 nm, Laserglow, Inc.) 
on an Olympus 1X7 1 microscope equipped with an Andor 
IXON (860) CCD camera. A 100-ms integration time was 
used for telomere DNA GQ and polyT measurements, 
and a 33-ms integration time was used for hairpin meas- 
urement. The stretching force applied to individual DNA 
tethers was altered by changing the linear position of a 
pair of rare-earth magnets mounted on a computer- 
controlled translation stage above the sample. The DNA 
extension and applied stretching force were determined by 
visualizing the magnetic bead with a blue LED (Thorlabs) 
mounted above the magnet assembly as described previ- 
ously (28). Before collecting smFRET data, the length of 
each DNA tether was monitored in real time as a function 
of applied twist to ensure that there was only a single 
DNA molecule attached between the bead and the 
surface. Molecules for which the extension changed as a 
function of twist represented DNA braids and were not 
used in the measurements. 

Raw movie files were first analyzed using in-house- 
written software (Labview and MATLAB), and then fit 
with the HaMMy (27) software to generate idealized 
FRET traces used for dwell time analysis. Single-molecule 
FRET histograms generated from combined FRET mag- 
netic tweezers data were compiled from all molecules 
analyzed at a particular force set point. Traces were first 
normalized for photo bleaching and then cropped before 
the photo bleaching event for HaMMy analysis. 
We focused our analysis on the equilibrium between the 
predominant mid-FRET (folded) and the low-FRET 
(unfolded) state. To this end, for our analysis, we 
used the idealized FRET traces generated by the 
HaMMy (27) fitting to remove the minority of time 
spent in the alternative high-FRET (folded) state. For 
dwell time analysis, idealized traces produced by 
fitting the smFRET trajectories using the HaMMy 
software package were used to calculate the average 
lifetime of the folded or unfolded states at each of the 
indicated forces. 

RESULTS 

Single-molecule FRET analysis of telomere DNA GQ 
structure in the absence of force 

The experimental geometry used in the present study of 
the force-dependent folding/unfolding of telomere DNA 
GQ structure includes attachment of duplex DNA handles 
to each end of a single GQ forming sequence (Figure 1A, 
right). A similar design has been used in recent optical 
trapping experiments (24-26), but it has not been 
analyzed using smFRET; therefore, we first set out to 
characterize the structure and dynamics of our model 
human single-stranded telomere DNA construct [Tel24, 
(TTAGGG) 4 ] embedded within a duplex DNA molecule. 
For these experiments, the Tel24 sequence was flanked by 
non-telomeric DNA extensions, which were hybridized to 
their respective complementary DNA strands harboring 
either a FRET donor (Cy3) or acceptor (Cy5) dye 
(Figure 1A, right). To minimize the possibility of the adja- 
cent duplex altering the folding properties of the GQ 



forming sequence, abasic sites were introduced between 
the Tel24 and the non-telomere sequence. Experiments 
investigating the structure of Tel24 in the absence of 
force were conducted on a prism-type total internal reflec- 
tion fluorescence (TIRF) microscope (29). Throughout 
this study, we define FRET = IaKJa + Id), where I A and 
Id are the background-corrected intensities of the acceptor 
and donor dyes, respectively. 

We first measured the folding properties of Tel24 over a 
wide range of NaCl concentrations (Supplementary 
Figure SI). As the NaCl concentration was gradually 
increased to 100 mM, we observed the emergence of a pre- 
dominant FRET population centered at FRET = 0.54, a 
second population centered at FRET = 0.33 and a minor 
population centered at FRET = 0.70 (Figure IB). 
Analysis of single-molecule trajectories revealed the 
Tel24 construct is in a dynamic equilibrium between 
each of these three FRET states (Figure 1C). Conversion 
of the predominant mid-FRET state to the high-FRET 
state required formation of the transient low-FRET 
state, as has been previously reported (21). Single-mol- 
ecule trajectories were fit to a hidden Markov model, 
which yielded idealized FRET trajectories (Figure 1C) 
(27). The idealized traces were used to generate transition 
density plots, which clearly indicated the majority of 
FRET transitions (70%, n = 2426/3475 transitions) 
occurred between the mid-FRET and the low-FRET 
states (Figure ID). The ratio of time spent in the 
unfolded (FRET = 0.33) and the predominant folded 
(FRET = 0.54) states yielded a K eq = 0.25, corresponding 
to a AG un foid = 0.81 kcal mol -1 , a value that is in good 
agreement with previously reported calorimetric studies of 
Tel24 unfolding in NaCl (30). Dwell time analysis yielded 
rate constants of 0.87 s _1 and 0.29 s" 1 for the folding and 
unfolding reactions, respectively, giving a K eq = 0.33 in 
reasonable agreement with the value obtained by fitting 
the smFRET histograms. Nuclear magnetic resonance and 
CD spectroscopy experiments on telomere DNA sub- 
strates in the presence of NaCl have demonstrated the 
strong preference of single-stranded telomere DNA to 
form the anti-parrallel GQ conformation under this 
folding condition (14,31). Thus, while smFRET data 
alone are not sufficient to conclusively determine which 
GQ structure is present, it is likely the predominant 
FRET = 0.54 conformation observed in the presence of 
100 mM Na + represents the anti-parallel conformation 
and the FRET = 0.33 state represents the unfolded 
conformation. 

In contrast to the Na + folding condition, the presence of 
KC1 produced a considerably broader distribution of 
FRET states, with increasing concentrations of K + 
stabilizing several distinct high-FRET GQ conformations 
(Supplementary Figure S2). The smFRET distribution for 
Tel24 folded in 100 mM KC1 is well fit by two Gaussian 
functions, with populations centered at FRET = 0.61 and 
FRET = 0.76. Consistent with this result, analysis of 
single-molecule trajectories in the presence of lOOmM 
KC1 revealed two distinct long-lived FRET states, which 
occasionally inter-converted via the obligatory low-FRET 
intermediate (Supplementary Figure S3). In addition, the 
presence of Mg 2 in the K + folding reaction appeared to 
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Figure 1. Single-molecule FRET analysis of Na -induced telomere DNA G-quadruplex folding in the absence of force. (A) Top left: Diagram of 
H-bonding network within a single G-quartet with a monovalent metal ion coordinated at its center. Bottom left: Schematic illustration of the 
anti-parallel G-quadruplex folding topology, with each planar G-quartet represented as a quartet of blue rectangles. Right: Schematic diagram of 
smFRET experimental setup [a model telomere DNA sequence ((TTAGGG) 4 , Tel24]. G-quartets highlighted in blue are flanked by two duplex DNA 
handles harboring a FRET donor (Cy3) and acceptor (Cy5) dye, respectively. Tel24 molecules were surface immobilized via a biotin-streptavidin 
linkage and imaged using prism-type total internal reflection microscopy. The antiparallel G-quadruplex structure solved by NMR is shown in the 
illustration (PDB code 143D). (B) Single-molecule FRET histogram derived from data collected on > 100 molecules is fit with three Gaussian 
functions centered at FRET = 0.33, 0.54 and 0.70. (C) A representative smFRET trajectory for Tel24 in the presence of 100 mM NaCl. Individual 
donor (green) and acceptor (red) dye intensities (top) and calculated FRET ratios (bottom) are plotted as a function of time. FRET trajectories were 
fit with a hidden Markov model yielding idealized FRET traces (magenta). (D) Idealized FRET traces from HaMMy fitting were used to generate a 
transition density plot indicating the frequency of each observed transition. Color bar on right indicates the number of transitions. 



have a slightly stabilizing effect on GQ structure, but did 
not significantly alter the smFRET histograms in the 
presence of lOOmM K + (Supplementary Figure S4). 
Taken together, our smFRET results are in good agree- 
ment with previous findings that demonstrated Na + ions 
promote homogeneous folding of telomere DNA GQs, 
while K + ions induced a more complex distribution of 
FRET states that are more thermodynamically stable 
than the Na + -induced fold. Based on these results, we 
conclude that the presence of duplex DNA flanking the 
Tel24 sequence does not significantly alter telomere DNA 
GQ folding and stability. 

Integrated fluorescence and magnetic tweezers 
spectroscopy 

Previously reported force spectroscopy experiments using 
optical traps achieved sub-nanometer spatial resolution by 
applying relatively large stretching forces (>20pN) to the 
system of study (32,33). Because many biologically im- 
portant structural transitions in nucleic acids and 
proteins are induced by much smaller forces, we have de- 
veloped a simple method for measuring sub-nanometer 



structural rearrangements within individual DNA mol- 
ecules held under a wide range of stretching forces (0.1- 
50 pN). In our integrated fluorescence and magnetic 
tweezers system, a single Tel24 sequence is embedded 
within a duplex DNA molecule specifically attached 
between a microscope slide and a magnetic bead 
(Figure 2A). The length of the DNA handle attached to 
the microscope slide was 29 bp (MOnm), while the length 
of the DNA handle attached to the bead was made con- 
siderably longer (15.7 kb, 5.4-um contour length) to 
prevent background signal introduced by the large 
magnetic bead entering the evanescent field produced 
near the surface of the glass slide by total internal reflec- 
tion of the excitation laser. The DNA construct is further 
modified with FRET donor (Cy3) and acceptor (Cy5) dyes 
on the DNA handles flanking the Tel24 sequence, as well 
as abasic sites between the duplex DNA handles and the 
Tel24 sequence. Stretching force was applied to the 
tethered DNA molecule with a pair of permanent 
rare-earth magnets mounted on a computer-controlled 
translation stage. Integrated fluorescence and magnetic 
tweezers experiments used objective-type TIRF, which is 
readily compatible with the presence of the magnet 
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Figure 2. Force biases the telomere DNA G-quadruplex folding/unfolding equilibrium. (A) Schematic diagram of the experimental setup for 
integrated fluorescence and magnetic tweezers measurements. A FRET-labeled Tel24 molecule is flanked by a short biotinylated duplex DNA 
handle on one side, and a second longer duplex DNA handle with terminal digoxigenin modifications on the other side. Individual Tel24 constructs 
were attached between a streptavidin-coated microscope slide and an anti-digoxigenin-coated magnetic bead in a magnetic tweezers setup and imaged 
using objective-type total internal reflection microscopy. Variable stretching force was applied to individual Tel24 molecules by translating a pair of 
rare-earth magnets held above the sample chamber. (B) Representative smFRET trajectories of a single Tel24 molecule held under three different 
forces (0.7, 2.5 and 5.1 pN). Histograms of the FRET values obtained at each force set point are shown on the right side of the FRET trajectories. 



assembly held above the sample chamber. We note that 
owing to slight differences in the optical components in 
the microscopes, smFRET values made by objective-type 
TIRF vary slightly from those measured using the 
prism-type TIRF setup described in Figure 1 [ie. mid- 
FRET (folded state) is 0.64 in low-force objective-type 
TIRF experiments versus 0.54 in prism-type TIRF micro- 
scope under zero force]; however, this variation does not 
impact our assignment of the mid-FRET and low-FRET 
states to the folded and unfolded telomere DNA GQ con- 
formations, respectively. 

In Figure 2B, representative smFRET trajectories are 
plotted at three different stretching forces for a single 
telomere DNA GQ in the presence of 100 mM NaCl. 
It is evident that at low forces (0.7 pN), the molecule 
spends the majority of the time in the high-FRET 
(folded) state consistent with our zero force measurements 
described earlier in the text. In contrast, a higher force 
(5.1 pN) substantially shifts the telomere DNA GQ 
folding equilibrium toward the low-FRET (unfolded) 
state. Notably, the application of small degrees of 
tension appeared to substantially reduce the number of 
transitions into the minority high-FRET (FRET = 0.70) 
conformation when compared with experiments per- 
formed in the absence of force; thus, we focused our 
analysis on the structural equilibrium between the 
low-FRET (unfolded) and predominant mid-FRET 
(folded) conformation. These data qualitatively demon- 
strate the capacity of our technique to detect nanometer- 
scale structural transitions within single telomere DNA 
GQs at low forces. However, to determine whether this 



approach can be used to extract quantitative information 
about force-induced structural transitions in DNA, we 
next performed control experiments on a model DNA 
hairpin construct (Supplementary Figure S5). Data 
taken on the model DNA hairpin were analyzed as a 
simple two-state system separated by a single energetic 
barrier, yielding an unfolding force Fi /2 — 12.1 pN 
(defined as the force at which the hairpin folded and 
unfolded states are equally populated) and a distance 
between the folded and unfolded state Ax = 16.2 nm. 
Furthermore, analysis of the force-dependent rate con- 
stants for folding and unfolding placed the transition 
state barrier 8.9 nm from the folded state and 6.1 nm 
from the unfolded state. These results are in good quan- 
titative agreement with previously reported optical 
trapping measurements on the same model DNA hairpin 
(2), demonstrating the ability of our integrated fluores- 
cence and magnetic tweezers system to extract quantitative 
information about DNA folding processes. 

Force dependence of telomere DNA GQ folding and 
unfolding 

We next characterized in greater detail the force depend- 
ence of the telomere DNA GQ folding/unfolding equilib- 
rium. For these experiments, we performed force 
titrations for 75 individual telomere DNA GQ molecules 
(for each force, a minimum of five different molecules 
were analyzed), and compiled the results into FRET histo- 
grams (Figure 3A, see Supplementary Figure S6 for 
complete data set). Analysis of the amount of time spent 
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Figure 3. Telomere DNA G-quadruplex folding equilibrium as a 
function of applied force. (A) Compiled single-molecule FRET histo- 
grams from data taken on 75 individual Tel24 molecules. For each 
force, a minimum of five different molecules were analyzed. Data for 
five representative forces ranging from 0.7 to 5.1 pN are shown (see 
Supplementary Figure S6 for complete data set). (B) Plot of the 
ln(K eq ) versus force. For each force set point at which folding/unfold- 
ing transitions were detected, smFRET histograms were fit with two 
Gaussian functions, and Keq for unfolding was calculated as the area 
under the low-FRET (unfolded) state divided by the area under the 
high-FRET (folded) state. Plot was fit with the expression 
ln(K eq ) = FAx/k B T — AG°/k B T, yielding an unfolding distance 
Ax = 2.7nm and a F 1/2 = 2.5pN. T-axis error bars are estimates of 
the standard error of the mean determined by bootstrapping analysis 
of smFRET data, and .Y-axis error bars reflect an ~10% systematic 
error in determination of stretching force. 



in the mid-FRET (folded) state versus the low-FRET 
(unfolded) state as a function of force provided a direct 
measurement of the effect of stretching force on the 
telomere DNA GQ folding equilibrium. We fit each 
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Figure 4. Force response of Tel24 low-FRET (unfolded) state and 
polyT DNA. The average FRET values of the Tel24 low-FRET 
(unfolded) state (open circles) and a 24-nucleotide polyT DNA con- 
struct (closed circles) are plotted as a function of force. 7-axis error 
bars are the standard deviation of the average FRET value measured 
for at least five individual molecules at each force set point, and x-axis 
error bars reflect an ~10% systematic error in the calibration of the 
applied stretching force. 



histogram with two Gaussian functions and calculated 
K eq for the unfolding reaction as the area under the 
low-FRET peak divided by the area under the 
mid-FRET peak. From the slope of the plot of ln(Keq) 
as a function of stretching force, we determined the un- 

2.5 pN 



1/2 



folding distance Ax = 2.7nm and 
(Figure 3B). 

We observed that the centers of the Tel24 low-FRET 
(unfolded) and mid-FRET (folded) populations shifted 
toward lower values as the applied stretching force was 
increased (Figure 3A). The two different Tel24 conform- 
ations fall within significantly different regions of the 
FRET response curve; thus, the magnitude of the 
force-induced A FRET for the low-FRET population 
reflects a substantially larger change in distance between 
the FRET probes than the AFRET of the mid-FRET 
population (Supplementary Figure S7). The low-FRET 
Tel24 conformation, previously identified as an obligatory 
folding intermediate during inter-conversion between 
distinct topological isomers (21), was suggested to reflect 
an unfolded DNA conformation that would be expected 
to respond to stretching force by gradually increasing its 
end-to-end distance, as observed in our experiments. To 
further explore the structural properties of this obligatory 
unfolded intermediate, we compared the force response of 
the low-FRET (unfolded) telomere DNA GQ conform- 
ation with a 24-nucleotide polyT DNA strand, which 
has been used as a model for unstructured single-stranded 
DNA in previous smFRET experiments (34). We find that 
under low stretching force conditions, the unfolded con- 
formation of the Tel24 construct is significantly more 
compact (higher FRET) than the polyT construct 
(Figure 4). This result demonstrates the G-rich Tel24 
unfolded state experiences an attractive force, perhaps 
derived from increased stacking interactions between 
purine bases or transient H-bonding, which promotes 
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telomere DNA GQ folding. Interestingly, the elastic 
properties of the unfolded Tel24 and polyT constructs 
converge within the force range (~8pN) that completely 
inhibits telomere DNA GQ folding within the 100-ms time 
resolution of these measurements (Supplementary 
Figure S6). We note that the gradual increase in extension 
of the unfolded Tel24 state with increasing force is not 
consistent with the presence of any long-lived structured 
GQ folding intermediates such as a DNA hairpin or 
triplex structure. These structured intermediates have 
been previously suggested to occur during telomere 
DNA GQ folding (35-37), but may be too transient to 
detect in our experiments. 

The force dependence of the rate constants for telomere 
DNA GQ folding and unfolding can also be analyzed 
from our single-molecule trajectories. Making the sim- 
plifying assumptions that along the reaction coordinate 
defined by the DNA-stretching axis, the position of the 
energy barrier is independent of the stretching force and 
the positions of the folded and unfolded energy wells are 
not substantially altered by the applied force, we modeled 
the lifetime of the folded state, ifoided. as given by 
tfoidedCF) = T folded ,o x exp(FAx :i: fold ed/k B T), where T fo i d ed,o 
is the folded lifetime at zero force, Ax^^ded is the 
distance from the folded to the transition state, k B is 
Boltzman's constant and T is temperature. Similarly, the 
lifetime of the unfolded state as a function of force can be 
expressed as T un(o[ded (F) = T unfo ided,o x expCFAxVfoided/ 
k B T). Fitting these expressions to a plot of the average 
lifetimes of the folded or unfolded states as a function of 
force yielded Ax^^ded = 2.1 nm, Ax^ed = 0.6nm, 
Tfoided.o = 4.1 s and T fo ided,o = 0.76 s (Figure 5A). The 
results of our kinetic analysis are in good agreement 
with the Ax value derived from the equilibrium data 
(Ax = Ax'" folded + Ax _t unfolded = 2.7 nm). In addition, the 
extrapolated zero force rate constants of kf ohl = 1.3 s -1 
and k un f oU i = 0.24 s -1 are in good agreement with our 
smFRET measurements in the absence of force (see 
earlier in the text). The short distance (0.6 nm) from 
the folded state to the transition state barrier suggests 
the telomere DNA GQ structure is stabilized by 



long-range contacts, and that disruption of just a few of 
the terminal base pairing contacts is sufficient to promote 
the unfolding of GQ structure (Figure 5B). 



DISCUSSION 

Here, we describe an integrated fluorescence and magnetic 
tweezers technique that combines objective-type total 
internal reflection microscopy with a 'classical' magnetic 
tweezers setup (28,38). Although several groups have 
described combined fluorescence-force instruments using 
optical traps to apply stretching forces (6,9), these tools 
remain inaccessible to non-specialists, owing to the tech- 
nical challenges associated with aligning and calibrating 
these advanced optical trapping systems. Furthermore, 
the integration of a high-power trapping laser with 
single-molecule fluorescence detection has the potential 
to negatively impact the photo physical stability of the 
fluorescence probes, thereby limiting the duration of the 
FRET measurement. By comparison, the technique 
described in the present study is a simple instrument to 
construct and maintain, and should therefore be a gener- 
ally accessible tool for a wide range of biochemists and 
structural biologists. Moreover, the use of permanent 
magnets to apply stretching forces obviates concerns of 
photo damage induced by a high-power trapping laser 
and provides a means to stably apply low forces for arbi- 
trarily long periods. Indeed, smFRET trajectories 
acquired using our approach routinely last for tens of 
minutes, facilitating measurement of a single molecule 
under a variety of experimental conditions. An additional 
advantage of hybrid force-fluorescence methods using 
magnetic tweezers is the ability to monitor torque-induced 
structural transitions by smFRET, as has recently been 
reported in studies of the torque-induced B-DNA to 
Z-DNA transition (7). We expect the capability of 
integrated fluorescence and magnetic tweezers systems to 
detect nanometer-scale structural rearrangements under 
conditions of low tensions and torques will be broadly 
applicable to a wide range of macromolecular folding 
studies. 




Force, pN Distance, nm 

Figure 5. Force-dependent rate constants of telomere DNA G-quadruplex folding and unfolding. (A) Average lifetimes (t) of the folded (high 
FRET) and unfolded (low FRET) states are plotted as a function of force. Plots are fit with exponential functions described in the main text, yielding 
As'unfoiijai = 2.1 nm and Ax*f 0 i<ied = 0.6 nm. (B) A schematic folding energy landscape for Tel24 in the presence of lOOmM NaCl. Cartoon model for 
the folded, transition-state and unfolded structures are superimposed on the folding energy landscape. The short distance between the transition-state 
barrier and the folded conformation suggests just a few base pairs are disrupted in the transition state structure. 
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In the present work, we have analyzed the force- 
dependent folding of a model human telomere DNA GQ 
in the presence of Na + . We find that the dynamic equilib- 
rium between the Na + -induced GQ folded and unfolded 
states is extremely sensitive to applied stretching forces 
between ~1 and 8 pN. By analyzing the effect of stretching 
force on the rate constants for telomere DNA GQ folding 
and unfolding, we have directly measured the position 
of the transition state barrier for GQ unfolding along 
the well-defined DNA-stretching reaction coordinate. 
Interestingly, we find the distance between the telomere 
DNA GQ folded state and the transition state barrier 
for unfolding is short (~0.6nm). Our finding is consistent 
with a recent report that used dynamic force spectroscopy 
(DFS) to reconstruct the distribution of rupture forces for 
a DNA GQ formed by a sequence upstream of the human 
insulin gene promoter (26). Fitting the DFS results to the- 
oretical models that predict the distribution of rupture 
forces assuming a single transition state barrier suggested 
the barrier for folding is within ~1 nm of the GQ folded 
state, in close accord with our present results. However, as 
with previous optical trapping experiments, which 
analyzed the mechanical stability of GQ DNA (24,25), 
the DFS measurements did not directly observe a revers- 
ible folding/unfolding equilibrium in the low force 
regimen as reported here. At present, it is unclear why 
the low-force 'hopping' behavior we report for the 
telomere DNA GQ in the present study is not observed 
in the optical trapping experiments, but the discrepancy 
may be due to differences in the sequences being used or 
challenges associated with detecting small extension 
changes at low stretching forces in an optical tweezers 
instrument. 

Force spectroscopy experiments with RNA and DNA 
hairpins revealed the transition state barrier for unfolding 
lies closer to the unfolded state than the folded state (1,2). 
This long distance to the transition state barrier for un- 
folding is expected for compliant structures the folding 
stability of which is derived primarily from local inter- 
actions (i.e. zipping of adjacent base pairs within a 
DNA hairpin). However, a short distance to the transition 
state barrier for unfolding, as measured here for the 
telomere DNA GQ, is indicative of a more brittle struc- 
ture the stability of which is reliant on critical contribu- 
tions of long-range interactions. This result can be 
explained within the framework of the nuclear magnetic 
resonance structure for the Na + -induced anti-parallel GQ 
fold (14), wherein guanine bases separated by as many as 
20 nucleotides must be brought together within a planar 
G-quartet configuration. Thus, disruption of just a few 
H-bonding interactions by mechanically pulling on the 
terminal bases of the telomere DNA GQ structure is suf- 
ficient to effectively destabilize the entire fold (Figure 5B), 
a finding that has direct implications on the molecular 
mechanisms of telomere-associated proteins and 
enzymes, which must resolve DNA GQs during telomere 
maintenance. Telomere DNA GQs are typically classified 
as DNA secondary structure; however, our studies reveal 
these unique DNA folds share important structural 
properties with RNA and protein tertiary structures, 
which have also been characterized by sub-nanometer 



distances to the transition state barrier for unfolding 
(1,39-41). Finally, understanding the critical contribution 
of long-range interactions in promoting telomere GQ 
folding stability should aid efforts in designing small-mol- 
ecule drugs that target and stabilize GQ structure to 
disrupt telomere homeostasis. 
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